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Abstract 



The Casimir energy of a non-uniform string built up from two pieces with 
different speed of sound is calculated. A standard procedure of subtracting 
the energy of an infinite uniform string is applied, the subtraction being in- 
terpreted as the renormalization of the string tension. It is shown that in the 
case of a homogeneous string this method is completely equivalent to the zeta 
renormalization. 
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I. INTRODUCTION 



In the standard setting of the Casimir effect problem, a bounded configuration space of 
the quantum field system is considered . However unbounded configuration space can be 
nonhomogeneous, i.e., it can consist from separate regions endowed with different physical 
characteristics. In the case of electromagnetic field such caracteristics are permittivity e and 
permiability /i of the media. It is obvious that the vacuum energy should depend on the 
configuration of these nonhomogenities. Hence the problem on finding out this dependence, 
i.e. on calculation of the Casimir energy, can be set. 

When calculating the vacuum energy it proves to be essential whether the velocity of 
light (in the general case the velocity of relevant quanta) is continuous when crossing the 
interface between regions with different characteristics. From a mathematical stand point 
the constant velocity of light implies that the coefficients of the higher (usually the second) 
derivatives in the corresponding dynamical equation are continuous functions in space. If it 
is the case then the calculation of the Casimir energy (or the Casimir forces) turns out to 
be practically the same as for an empty space with a perfectly conducting shell having the 
shape of the interface between media 0-0. 

The discontinuity of the speed of light on the boundaries between different regions results 
in considerable mathematical difficulties. Certain results have been obtained here only 
recently. The Casimir energy was calculated for a compact dilute dielectric ball proceeding 
from the Green's function formalism and employing the "naive" zeta function technique for 
removal of the divergences without detailed procedure of analytic continuation |9|-[T2|. It 
was assumed that the difference of the light velocities inside (ci) and outside (02) is small 
and only the term proportional to {c\ — C2) 2 in the Casimir energy has been kept. In Ref. [13 
it was shown that the heat kernel coefficient 02, responsible for the pole contribution to the 
Casimir energy, vanishes in this approximation (it proves to be proportional to (ci — C2) 3 ). 
This implies that the complete zeta regularization should provide a finite answer in this 
problem. 

In calculations accomplished in Refs. [§,1C,1^] a key role was played by the so called 
contact terms which cannot be incorporated in the standard zeta function technique at least 
in a straightforward way. Their physical origin is still unclear. 

Thus it is worth investigating these problems in the framework of a more simple model. 
Such a model, preserving all essential features of field theory, is a string built up from two 
pieces which have different velocities of sound (the velocities of vibration propagation). In 



papers by I. Brevik and co-authorths [14-Tq] a piecewise uniform string, with the same 



velocity of sound at each of its constituents, has been studied. Such string model is a 
one-dimensional analog of the electromagnetic field in media obeying the condition efi = 
c~ 2 , where c is the constant velocity of light. In this model the Casimir energy has been 
investigated in detail both at zero and finite temperature. 

The present paper seeks to calculate the Casimir energy of a nonuniform open string 
of total length R, which is built up of the pieces of length r and R — r having the sound 
velocities vi and v-i respectively. At the joint point (x = r) the continuity of the string shape 
and its tangents is imposed. The string ends x = and x = R are subjected to the Dirichlet 
or Neumann boundary conditions. The equations defining the eigenfrequencies of the string 
are derived. This material is presented in Sect. II. Proceeding from this in Sect. Ill the 
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Casimir energy of the string is calculated by making use of the mode-by-mode summation 
method which relies on the contour integration. The renormalization of the string 

tension by subtracting the energy of an infinite uniform string enables one to remove the 
divergences in a unique way. It is shown that in the limit R — > oo the Casimir energy of the 
string does not depend on the boundary conditions imposed at x — R and is expressed in 
terms of the polylogarithm function. In the Conclusion (Sect. IV) the results obtained are 
discussed briefly. In the Appendix it is shown that the subtraction procedure employed in 
the paper is equivalent to the zeta renormalization when considering a usual uniform string. 



II. CLASSICAL DYNAMICS OF A NON-UNIFORM STRING 



The model we are considering is a string built up from two pieces of length r and R — r 
and endowed with the sound velocities V\ and t>2, respectively. The string is stretched along 
the X 1 axis (which we shall denote by x in what follows) and can vibrate in D — 2 transverse 
dimensions X 2 . . . X D . The both pieces of the string satisfy the wave equation 



1 d 2 X i d 2 X l 



v] dt 2 



dx 2 



0. 



J 



1,2, 



2, 3,. ..,D. 



(2.1) 



At the ends the string is subjected to the boundary conditions, which we shall take to be of 
either Dirichlet or Neumann type. Consequently, we shall consider four cases: 
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(2.2) 



At the joining point x = r we impose the continuity conditions 

lim X%x) = lim X%x), 



lim 

x— »r_ 
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dx 



lim 
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'dX i 
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(2.3) 



The wave equations (|2~T| ) together with the boundary conditions ( |2.2| ) and the continuity 
conditions ( p. 3D lead to the equations which determine the eigenfrequencies of the string u 

fx(u) = sin(acj) + 5v sin(/5cj) = 0, (DD), 

f 2 (u) = cos(aco) - 5v cos(/^) = 0, (DN), 

f 3 (w) = cos(au) + 5v cos((3uj) = 0, (ND), 

U(lu) = sm(acu) - 5v sm(pto) = 0, (NN), (2.4) 
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where 

. v\ — v 2 r R — r r R 

ov = , a = 1 , and p — — 



Vi + v 2 Vi v 2 Vi v 2 

Let us note, that a > and a > (3. 

III. GROUND STATE ENERGY OF THE STRING AND ITS 

RENORMALIZATION 

When calculating the Casimir energy, the main problem encountered here is the removal 
of the divergences, which arise inevitably in such calculations. In quantum field theory fI7 



this procedure is accomplished in the course of the renormalization of the physical param- 
eters, which specify the field theory model under investigation (the masses of the quanta, 
coupling constants an so on). In order for this procedure to be correct from mathematical 
stand point, divergent expressions should be regularized, the regularization being removed 
after renormalization. 

We define the renormalized Casimir energy in a standard way [[l 

E c = \ ( w « _ m n) , (3.1) 

where the u's are the eigenfrequencies of the string, i.e., the roots of the frequency equa- 
tions (|2.2|), and the cJ's are the roots of the same equations but in the limit when the total 
string length R and r tend to infinity, and the difference of the velocities 8v tends to zero. 
Hence in this limit the initial string becomes a uniform infinite string. In quantum field 
theory this implies the removal of the Minkowski space contribution |l] . In the Appendix it 
is shown that the method of calculation of the vacuum energy, outlined above, is completely 
equivalent to the zeta regularization in the case of usual uniform string. 

It is convenient to represent the sum (|3.1| ) in terms of the contour integral 



where the function f(z) is any function fi(z), I = 1, ... ,4 defining the frequency equations 



in (2.4), and / is obtained from / when passing to the limit described above. The contour 
C encloses all the positive roots of the equation f(u) = 0. This contour can be deformed 
into the semicircle C\ of radius A in the right half-plane and the interval of the imaginary 
axes (— iA,iA). When the radius A is fixed, the contour C encloses a finite number of the 
roots of the equation f(z) = 0. The sum of these roots is obviously finite. Hence the radius 
A is a regularization parameter, and taking the limit A — > oo means the removal of the 
regularization. 

The subtraction in Eq. (|3.2j ) accomplished with the function f(z) can be interpreted 
as the renormalization of the parameters which specify the classical energy of the string. 
From the general consideration one can assume that this energy is proportional to the total 
length of the string R, i.e., E d = T ■ R, where T is a dimensional parameter [T] = L~ 2 . It 
is determined by a concrete string model that should certainly involve nonlinearities. For 
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example, in the relativistic string model |19j this parameter is the string tension. Due to 
the quantum corrections the string tension should be renormalized in the following way 



rr cn = T + £2 1 r dz d^ 
R Am Jc dz 



After removing the regularization (A — > oo) the integral in ( p.3| ) obviously diverges. However 
it is assumed that this divergence is canceled by the respective infinity of the "bare" string 
tension T. As a result the renormalized (physical) string tension T ren proves to be finite. Here 
we follow the standard procedure of removing the divergences in quantum field theory fl7 |. In 
view of the oscillating behavior of the functions fi{z), I = 1, . . . , A on the semicircle C\ their 
asymptotics for R — > oo are simply these functions themselves, i.e., fi(z) = fi(z), z e CV- 
Therefore the integration along this part of the contour C does not give contribution into 
the Casimir energy ( |3.2| ). On the imaginary axes the asymptotics of the functions fi(z) when 
R, r — > oo can be find easy 

m^—, «=-+—. (3.4) 

2 vi v 2 

As a result we obtain the final expression for the Casimir energy 

n — 9 r°° 

E c = ^— dylnMy)}, Z = l,...,4 (3.5) 
2n Jo 

with the functions hi(y) given by 

h x {y) = 1 -e~ 2ay + 25v e- ay sinh{(3y), (DD), 

h 2 (y) = 1 + e~ 2ay - 25v e~ ay cosh([3y), (DN), 

fa (y) = l + e~ 2ay + 25v e~ ay cosh((3y) , (ND) , 

h A (y) = 1 - e~ 2ay - 25v e~ ay sinh(/?y), (NN). (3.6) 



As a simple consistency check, we apply the formulas (|3.5| ) and to the uniform 

string {yi = v 2 = v ) of length R with the expected result 



D-2 , , / 9 B,N n(D-2) 
E c = / dy\n(l-e- 2 - y ^ - 



2tt Jo " \ > 2AR 

in the DD and NN cases and 



£c = — /„ * ln ( 1+ 



A8R 



-v 



in the DN and ND cases. 

Let us consider the string configuration, when the second piece of the string becomes 
infinitely long (i.e. R tends to infinity, while r remains fixed). Then we have 

K(V) = h?(y) = 1-Sv e~ 2 ^ = h-(y), (DD), (DN), 

hf{y) = h?(y) = l + 5v e" 2 ^ = h+(y), (ND), (NN). (3.7) 



5 



These formulas imply that when R — > oo the Casimir energy becomes independent of the 
type of the boundary condition imposed at x = R. Certainly, this is an appealing property 
of the vacuum energy of the string under consideration. For R — > oo it is also possible to 
obtain the explicit form of this energy 

D - 2 r°° 
E ± = ^— dyln[h ± (y)} 

Z7T JO 

D — 2 f°° ( ara' 

= / ay In 1 ± ov e "i 

D_2 ^Li 2 (±^), (3.8) 



where 



47r r 



Li,(^) = E^7' \ z \ <:L » 
fc=i K 

f-i. 



r(z/) jo e* — ^ 



2 /-oo t v ~ dt 

> 0, |arg(l -2)| < 7r 

e" — z 

is the polylogarithm function ||21|| . 

Our consideration of non-uniform string can be generalized in a straightforward way to 



the finite temperature. It can be done by the substitution [22 



dy^2n9Y!s(y-n n )dy, (3.9) 

n=0 

where 6 is the temperature, Q n = 2im8 are the Matsubara frequencies, and the prime means 
that the term n = is taken with half weight and 



IV. CONCLUSION 

The non-uniform string with finite r and R — > oo is, in some sense, an analogue of 
the radial part of the problem concerning the calculation of the vacuum energy of a pure 
dielectric compact ball placed in a homogeneous unbounded medium or cavity in such a 
medium fHH. The main lesson of our consideration is the following. When defining the 



counter term in Eq. (3T), in order to renormalize the Casimir energy, we take not only the 
limit R, r — > oo but also put 5v — > 0. It means that all the physical parameters, specifying 
the problem under study, should be involved in determination of the counter term. 

It is worth noting that we do not use here the contact terms in order to get a finite result 
for the Casimir energy. 
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APPENDIX: RIEMANN ZETA FUNCTION IN THE UNIFORM STRING 

MODEL 



The representation of the spectral zeta functions in terms of contour integrals, with the 
integrands being the relevant frequency equations, is of wide use [|!3l,p3| This representation 



is a direct application of the principle of argument theorem from the complex analysis fL8| , |24 
However this theorem alone does not afford the required analytical continuation of the zeta 
function. The details of this technique can be clearly demonstrated by considering the zeta 
function for the usual uniform string. In this case we are dealing with the Riemann zeta 
function the analytic continuation of which is well known ITSLEl 



The eigenfrequencies of the uniform string of the length R with the Dirichlet boundary 
conditions 

mr 

lLl n = —, 71=1,2,... (Al) 

ri 

are the positive roots of the equation 

srn{u)R) = 0. (A2) 

For simplicity the velocity v is taken to be 1. Let us define the zeta function in the problem 
at hand by the standard formula 

oo 

COO = Ra>l. (A3) 

n=l 

When R = ir we obviously have the Riemann zeta function. 

Now we represent the sum (|A3|) in terms of the contour integral 



«'> = ^/ c ^3S h W* < A4) 



where the contour C encloses all the positive roots of Eq. (|A2[ ). The point z = is the 
branching point of the function z~ s = exp(— sin 2;) in Eq. (A3). Obviously this point should 



be left outside the closed contour C in order to integrate one valued function. For the 
logarithm we take as usual the branch acquiring real values on the positive half-axes < 
31 z < 00. Let us deform (still formally) the contour C in a semicircle of infinitely large 
radius laying in the right half-plane of the complex variable z and close it by imaginary axes 
—00 < Qz < 00. For 3?s > 1 the integration along the semicircle can be dropt. As a result 
one gets 

C(s) = - sin i^—j J o dyy- s —\n[smh(yR)] (A5) 

R /tts\ f°° 

sin — / dyy~ s coth(yR) . (A6) 



71 V 2 J Jo 



The integral in Eq. (AS) converges at the upper limit only for 3?s > 1 and at the lower 



limit for 3? s < 0. The latter is due to the singular behaviour ~ e s of the integral along 
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the semicircle of infinitely small radius e at the origin. This contribution vanishes only for 
5is < 0. Thus, in the general case the contour C in Eq. (A4) cannot pass through the point 
z = 0. 

However this contour can be retained if we take, instead of Eq. (|A2|) , a new frequency 
equation 

4ir = °- (A7) 

which has the same positive roots as Eq. (|A2|), but the point z = does not satisfy it. As 
a result we obtain, instead of Eq. 



« s ' = 7 sta (y)r* r *( coth(!/B) -^)' (A8) 

This integral representation of the zeta function is defined for 1 < dts < 2 and gives in fact 
the same sum (|A3|) but with explicit substitution of the roots ([AT]). Really, let us present 



the function coth(yR) in Eq. ( |Aq ) as a series (see, for example, Eq. 1.4.21.4 in f25| ) 

1 00 1 

cothz = - + 2zJ2- — . (A9) 

After that integration in ( [AS] ) can be done 

y l " s dy 



C( S ) = -sin - )2RY: 



R . 


/7TS\ 


= — sin 




7T 




R . 




= — sin 




7T 


\ 2 



2 



n=l 



2 J ^Jo y 2 R 2 + tt% 

r (f) r KHg(i) 

nit' 



E(i) • (MO) 



Thus the contour integration of the frequency equation ( |A2|) or (|A7|) gives us nothing new 
as compared with the sum ( |A3| ) with explicit frequencies (|Al|) . 

As known |l8| , p4]| analytic continuation of the series ( |A3|) with R = tt to the region 
9fts < 1 is provided by the formula 

Here CrI- 2 ) is the Riemann zeta function, and the path of integration encircles the real 
positive axes < 3ft z < oo. When 3? s > 1 the integration contour around the origin can be 
drawn to the point, and the remaining integral gives 



oo ^ s 1 



CR(a) = pfr/ -^-rdx. (A12) 
1 (z) Jo e x — 1 

Expanding the function l/(e x — 1) in powers of e x and integrating each term we arrive at 
the sum 
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OO 1 

E 



Cr=U- »s>1. (A13) 

If 3ft s < 1, the integral around the origin in Eq. ( All ) does not vanish, and it should be 
taken into account when reducing the contour integral representation ( All ) to the ordinary 
integral (see, for example, p4fl). 

The Riemann zeta function defined in the entire plane s, safe for the point s = 1, by the 
contour integral ( [All] ) obeys the Riemann reflection formula [|T8| , p4 



Cr(s) 



(27T) 



7T 



sin 



TTS 

y 



r(l-s)C R (l-s) 



This formula relates the values of the zeta function left and right from the line 3ft s 
the point s = 1 the function Cr( s ) nas a simple pole 



feOO 



+ 7 



(A14) 

: 1. At 

(A15) 



where 7 is the Euler constant. It is this pole that prevents the use of the definition (|A13|) 
for 3ft s < 1. 

Thus we see that even for analytic continuation of the simple sum ( |A13| ) a special choice 
of the path and the integrand is needed. Direct summation of the frequencies by contour 
integration does not afford analytic continuation. 

However if we define the zeta function for a uniform string by subtracting the vacuum 
energy of an infinite string from Eq. ( |A13| ), as it is usually done when calculating the Casimir 
energy (see Eq. (|3.1D ), then the resulting zeta function proves to be the analytic function 
for 3ft s < and exactly that function which follows from Eq. ( All ) or from the reflection 
formula ( |A14j) , that is the same. 

So let us define the zeta function in the following way 



COO = E ("C ~ 



(A16) 



where the frequencies uj are defined by Eq. (|A2|) when R —>■ 00. The corresponding limiting 
form of this equation can be found unambiguously upon rotating the integration path to the 
imaginary axes 



lim sinh(yi?) 



R—>oo 



Now we apply Eq. ( |A5| ) to the both terms in definition ( JA16 ) 



COO 



sin 



7T 



TTS 

Y 



-s-l 



In 1 



-2yR 



dy. 



(A17) 



(A18) 



Here the integration by parts has been done, that is legitimated for 3ft s < 0. Expanding the 
logarithm in ( |A18|) in power series and integrating each term we obtain for 3ft s < 
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w N (2vr) s /7TS\ „. , ~ 1 

c(s) = L_L sm (_) r(1 _ s) | ; _ 



(2 ^sin(^)r(l- S )C(l- S ). (A19) 



It is exactly the same result that follows from the Riemann reflection formula (|A14 ). 



Thus the analytic continuation of the Riemann zeta function (A13) into the region 3ft s < 



1, which is given by Eq. ( |A11| ), implies in essence that for 3ft s < 1 one considers the initial 
series ( |A13|) with subtraction ( |A16| ). Certainly for 3ft s > 1 the new definition of the zeta 
function ( |A16| ) cannot be reduced to the initial series ( |A13| ), hence it cannot play the role 
of Eq. ( |A11|) , i.e., it cannot afford, in a rigorous mathematical sense, analytic continuation 
of the series ( |A13| ) to the region 3ft s < 1. The new definition ( |A16| ) simply "guesses" the 
result of analytic continuation of the Rieman zeta function ( |A13| ) to this region. 

To our opinion, this consideration shows, in a simple and clear way, the relationship 
between the analytic continuation, providing the mathematical basis of the zeta renormal- 
ization technique, and the method of subtraction or counter terms, widely used by physicists. 

Summarizing we arrive at the conclusion that the calculation of the Casimir energy of a 
non-uniform string in the present paper is completely equivalent to the zeta regularization 
in the case of a usual uniform string. 
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